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ABSTRACT 

Massive black holes appear to be present in the nuclei of almost all galaxies, 
but their genesis and evolution are not well understood. As astrophysical black 
holes are completely characterized by their masses and spins, the observed joint 
distribution of these quantities contains important clues to their history. We 
examine the coevolution of mass and spin in binary merger growth scenarios. We 
find that holes are typically spun down by mergers. Rapid rotation results only if 
the binary's larger member already spins quickly and the merger with the smaller 
hole is consistently near prograde; or, if the binary's mass ratio approaches unity. 
If, as some observations have suggested, observed black holes spin rapidly, then 
this limits the importance of merger scenarios for the growth of black holes. 

Subject headings: black hole physics — gravitation — galaxies: active, nuclei — 
quasars: general 



1. Introduction 

Black holes span a wide spectrum of masses: the case for stellar mass holes (M ~ 10 Mq) 
in the field [e.g., Bailyn et al. (1998)] and supermassive holes (M ~ 10^ — IO^Mq) in galactic 
bulges [e.g., Ferrarese (2002); Kormendy & Gebhardt (2001)] is extremely strong; tantalizing 
evidence suggests middleweight holes (M ~ 10^ — 10*^ ^0) as well (Colbert & Mushotzky 
1999; Colbert & Ptak 2002; Gebhardt et al. 2000; van der Marel 2001). Stellar mass holes 
likely form in stellar collapse; the origins of more massive holes remains mysterious. Such 
holes could form in the collapse of massive gas accumulations; they could grow from smaller 
holes by accretion; they could grow by capturing stellar mass bodies; and they could grow by 
repeatedly merging with holes of comparable mass. Any or indeed all of these mechanisms 
could contribute to the growth of a given hole. 
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A black hole's spin may help identify which scenario most strongly impacted its recent 
history. Since spin likely drives outflows and jets in active galaxies, and since jets are 
presumed to align with black hole spin (Rees 1978), spin may provide an observational probe 
of a hole's recent growth (Merritt 2002). We examine how spin and mass coevolve in mergers. 
Binaries will form following galaxy mergers (Begelman, Blandford, & Rees 1980), and may 
harden to the point that gravitational-wave (GW) emission drives its members together. 
Eventually, they encounter the last stable orbit (LSO) , and then plunge and coalesce into a 
single hole. Our goal is to understand the mass and spin of this remnant hole. 

For nearly equal mass holes, this is extremely difficult: we must model the spacetime 
dynamics of the transition from a binary to a single black hole, accounting for both holes' 
spins and the radiated energy and angular momentum. A proper analysis requires mature 
numerical relativity codes [see, e.g., Lehner (2001)]. The problem is simpler for small mass 
ratio, q = 1112 /mi = m/M <S 1. This binary is well described as a test particle orbiting a 
black hole. GW emission shrinks the small hole's orbit to the LSO, whereupon it plunges 
into the large hole. Neglecting the final emission of radiation after the plunge, the hole 
evolves simply: its mass adds the small body's energy at the LSO, its spin adds the LSO 
angular momentum. 

Because we only need global "conserved" quantities, this description works surprisingly 
well even for rather large mass ratio. Post-Newtonian analyses (Blanchet 2002; Buonanno 
& Damour 1999; Damour 2001) show that finite mass ratio typically changes the LSO and 
its orbital constants by a factor of order 1] = mM / (m + M)^ = q/{l + g)^ < 0.25. The error 
due to the test particle description is < 0.3 for q < 0.5. We also may safely neglect the 
energy and angular momentum radiated in the final merger: although its GW luminosity 
may be large, its duration will be very short. The mass carried off in this phase, for example, 
is AM ~ (0.01 - 0.1)Mg2 (Davis, Ruffini, Press, & Price 1971; Sasaki & Nakamura 1982). 
Neglecting this radiation incurs an error that is less important than other errors built into 
our approximations, and rapidly becomes negligible for small mass ratio. Likewise, the small 
hole's spin can be neglected: since a hole's spin scales with its mass squared, spin will be 
less important than the orbital angular momentum, provided we exclude q > 0.5. 

We set the speed of light c and Newton's constant G to 1; a useful conversion factor 
is 1 A/q = 1.5 km. Our binary has masses M and m; the larger hole has mass M and spin 
|S| = aM = aMc/G] < a < M. Vectors are written in boldface; hatted quantities have 
been made dimensionless by dividing out powers of mass — e.g., d = a/M. 
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2. Orbital constants and black hole evolution 

On reaching the LSO, the smaller hole plunges into the large hole, carrying along its 
orbital constants — energy E, angular momentum parallel to the spin Lz, and "Carter con- 
stant" Q. The Carter constant separates the equations of motion in a Hamilton-Jacobi de- 
scription of black hole orbits [e.g., Misner, Thorne, and Wheeler (1973) (MTW)]. It is essen- 
tially just the "rest" of the orbit's angular momentum: to very high accuracy (Glampedakis, 
Hughes, Sz Kennefick 2002), one can describe the binary as having an angular momentum 
L"^ — Q + L^. We treat Q as L^, angular momentum projected into the eqiiatorial plane 
(perpendicular to the spin). This treatment resonates with the theory of orbits in axisym- 
metric potentials: Q is a relativistic analog of the "3rd integral" J3 [cf. Binney & Tremaine 
(1987)]. Treating Q as is exact for orbits of non-spinning holes; for maximal spin, the 
error is less than a few percent (Glampedakis, Hughes, & Kennefick 2002). E includes the 
rest mass of the orbiting body: bound orbits have E/m < 1, unbound orbits E/m > 1. 

The LSO is a one-parameter set of orbits: each orbit has radius r and constants 
{E, L^i Q) determined by the inclination angle t, defined as 

cos i = L2/ yjL\^Q = H. (1) 

This angle is very useful: detailed studies (Hughes 2001) show it remains practically constant 
during inspiral, so a distribution /(/i) describing an ensemble of binaries at formation likewise 
describes that ensemble at plunge. To find the constants for circular orbits at plunge, solve 

i? = 0, i?' = 0, i?" = 0; (2) 

prime denotes djdr^ and the "potential" R is given by 

R = {E{r'' + a^) - aL,\ ' - A(r) [r^ + (L, - aEf + Q\ , (3) 

where A(r) = — 2Mr + . This potential describes the orbit's radial motion; see MTW, 
Chap. 33. The LSO is bounded by prograde {\i = 1) and retrograde (/x = — 1) equatorial 
orbits, with constants (Bardeen, Press, & Teukolsky 1972) 

^Lso = tlso/M 

= 3 + ^2 T V(3 - ^i)(3 + Zi + 2Z2) , (4) 

1 - 2f 2 ± av^ ^ ^ 

Elso = £;Lso/m =-====== , (5) 

Vl — o?;^ ± 2mr 

f T / ,^ , 1 =F 2af ^ a^t;" 

Llso = WmM = ±r.-^===. (6) 
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where v = ^/l/f, and 

Z, = i+(i-a2)i/3[(i + a)V3 + (i_a)V3] , (7) 

- [3a2 + . (8) 

The upper sign is for prograde orbits, the lower for retrograde. We now drop the "LSO" 
subscript, since we always refer to these quantities at the LSO. 

Using the bounding initial guesses, it is straightforward to solve Eq. (2) using 

Newton's method (Press et al. 1992) for the constants. The results are surprisingly well fit 
by a simple rule: letting ^ stand for r, E, or L, 

~ \U\ + ^(/^ + l)(^pro - ■ (9) 

Using Eq. (9) for L, one builds L^(^) ~ L(/i)/i, L±{jj,) ^ L{jj,)y^l — jj,"^. These fits are 
extremely good for small spin, and induce errors of about 5 — 10% for a ~ 1. 

We should also consider eccentricity e. Major mergers are well described by a circular 
LSO (eccentricity rapidly bleeds away in the formation of a tight binary and by the inspiral), 
but minor mergers will have significant eccentricity. We find that a black hole growth history 
barely changes when eccentricity is taken into account: almost identical growth histories are 
obtained with e and with e = 1. We will confine our discussion to circular orbits. 

It is now simple to compute a remnant's properties. Before merger, the large hole has 
mass M and spin S = aM'^ along the z axis. After merger, the remnant has mass and spin 

M' = M[l + qE{d,ii)] , (10) 
S', - M2[a + g4(a,Ai)] , (11) 
S'^ = qM^L^{d,l^). (12) 

The remnant hole is inchned at an angle relative to the original hole, and has spin d': 

A0 = arccos (^S'J S'^^ + S'lj , (13) 

a' = ^JS',^ + S'l/M''' . (14) 

3. Results 
3.1. Single major merger 



We now examine the remnant's properties following a single merger, choosing q and 
computing {E, L^jQ) as functions of the larger hole's spin d and the inclination cosine /i. 



- 5 - 



We do not yet use the approximation (9), but instead solve Eq. (2) numerically. We then 
use Eqs. (10) - (14) to describe the remnant. 

Two examples, q — 0.1 and q — 0.5, are shown in Fig. 1. We show a' as a function 
of d and /i. Consider the q — 0.1 results first. For a broad range of d and /i, the remnant 
spins relatively slowly. In many cases, d' < d: over much of the parameter space, this hole is 
spun down by the merger. Rapid rotation follows only if the larger hole was already spinning 
rapidly and plunge occurred at shallow inclination (// — 1). Nearly nonspinning remnants 
form near jj, ~ —1, d ~ 0.4: the retrograde orbit cancels the hole's spin. This requires 

? - ?crit = a/4et(a) ^ 0.23 . (15) 

When q < q^rit, the spin orientation changes very little by the merger. By contrast, when 
q > Q'crit, the spin is overwhelmed by the plunging body, and the orientation aligns with the 
plunge angle. Whereas about half of the parameter space for q = 0.1 leads to remnants with 
spin d < 0.4, the half-area contour when q — 0.5 is at spin d ~ 0.8. 

The slow spin of most remnants spin is simply understood. Angular momentum at 
the LSO depends strongly on inclination — cf. the fit (9). The magnitude |L| is small for 
prograde orbits and large for retrograde orbits. Thus |L| is smallest when it tends to augment 
the spin (S and L nearly parallel) and is largest when it tends to cancel (S and L nearly 
antiparallel) . On average, the hole tends to spin down. This tendency breaks at large g, 
when L overwhelms S. The spin of the remnant is then dominated by the orbit at plunge. 



3.2. Repeated minor mergers 

After a hole grows above a certain mass, its subsequent spin evolution may be stochastic. 
This limit may be described by a Fokker-Planck equation, combining the secular and diffusive 
changes in the hole's characteristics. We treat In M and a as independent variables, so that 
the distribution function (per unit volume a space) /(a. In M) evolves via 

where di is the i-th component of a, 

and angle brackets mean to average over its distribution [see, e.g., Pathria (1972), Lifshitz 
& Pitaevski (1980)]. We restrict our attention to a population of black holes that are born 
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with a specific mass and spin. This solution can be used to integrate over a broad initial 
population. We also will limit our quantitative analysis to an isotropic distribution of (small) 
merging holes; generalization to an anisotropic distribution is straightforward though lengthy. 

We rewrite Eq. (16) in spherical coordinates attached to the initial spin direction, 
(a, 9, (f)). For the isotropic case, symmetry dictates that the only non-zero coefficients are 



R = 

Equation (16) then becomes 

df 
dlnM 



Aa 



AlnM 

(Aa)^ 
AlnM 



AlnM 



(18) 
(19) 



+ 



l_d_ 

0,2 da, 

2 sin 6 86^^^ 86 



ld_ 

d 



df 



(20) 



Evaluating the coefficients requires that we expand Eqs. (10) - (14) to leading order in q. 
For each merger, the changes in d and 9 satisfy 



5a — q Lz{d,iJ,) — 2d + 0{q'^) , 

69 = %L{d,ii)^l- ii^ + 0{q^) . 
a 



(21) 
(22) 



Combining this with 5\n.M = qE{d,ii) yields 

[L^(a,//) - 2a] 
E{d,fj,) 

'[4(a,/x) -2a]2 



R 



E{d,n) 



(q) / [L{d,fi)^/1 - fi 



E{d,fi) 




(23) 
(24) 
(25) 



(q) is the typical minor merger mass ratio. In the approximations, we put E — 1 and expand 
in powers of a; this is permissible for d < 0.9. The Fokker-Planck approach works well for 
{q) < 0.3, and is easily supplemented by direct sums (as in Sec. 3.1) for major mergers. 

The Fokker-Planck coefficients can be used to understand semi-quantitatively the evo- 
lution of holes that grow through minor mergers. The change in a is dominated by the 
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resistive term for R > (D|| + D±)/a, which is the case for a > 2(g)^/^. In this case, the 
fluctuations about the average evolution are small. This average, secular evolution can be 
found by using the definition of R and integrating Eq. (23) to find 

d{t) = a(0) [M{0)/M{t)f^ . (26) 

If the power were 2 rather than 2.4, this equation would tell us that the hole's original spin is 
preserved while the mass grows. The spin dies away somewhat faster because the magnitude 
of the change for retrograde captures is larger than that of prograde captures. 

The orientation evolves diffusively. After growing by AM, the typical misalignment is 



{Se) ~ y/2D^AlnM/a^ ~ 2.7 J— . (27) 

V M 

Significant realignment in a single merger occurs only if g/a ~ 0.3, in accord with Eq. (15). 



3.3. Rapidly spinning holes 

We have argued that a hole is unlikely to acquire rapid spin through mergers. Turn 
to the converse problem: if we believe that a black hole is spinning with a ^ 1, how many 
minor mergers can it have experienced? To be specific, suppose that the hole is born through 
gravitational collapse or is spun up by accretion to a ~ 1. If we assume that the hole spun 
down to its current value by capture from an isotropic distribution of minor mergers, then 
as da/ din M ~ —3 for a ~ 1, the mass acquired must satisfy 

AM < (1 - a)M/3 . (28) 

Rapidly spinning holes are extremely unlikely to have suffered a recent merger; their spin 
must either be original or due to accretion. 



4. Discussion 

We find that the remnant of a major merger is rarely rapidly rotating: rapid rotation 
follows only if the larger binary member spun rapidly before merger and the plunge was nearly 
prograde; or, if the binary's mass ratio 5 ~ 1. Given the variety of black hole masses seen 
in galaxies, mergers with 5 ~ 1 should be rare; given the small volume of parameter space 
leading to rapid rotation, serendipitous configurations leading to a rapid rotation should also 
be rare. Mounting evidence, mostly from spiral galaxies, suggests that in many cases massive 
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black holes nonetheless rapidly rotate [e.g., Wilms et al. (2001); Elvis, Risaliti, & Zamorani 

(2002)]. Our results strongly suggest that this spin cannot come from mergers [e.g. Wilson 
& Colbert (1995); Kauffmann & Haehnelt (2000)] but, instead, is consistent with the view 
that black hole mass (and spin) is assembled radiatively [e.g. Small & Blandford (1995); Yu 
& Tremaine (2002)]. 

The spin evolution of a hole that grows by repeated minor mergers is neatly described 
by a Fokker-Planck equation [Eq. (16)], taking a particularly simple form if the mergers 
arise from an isotropic cluster. In this case, the evolution has a secular component, which is 
approximately described by a "doctrine of original spin" [S = aM"^ remains roughly constant 
while M grows; cf. Eq. (26)], and a diffusive component, with a spectrum of fluctuations 
governed by the coefficient D|| [cf. Eq. (24)]. This is in accord with recent work on models 
to grow intermediate mass black holes in clusters (Miller 2002). 

Finally, we predict the typical angle (66) by which a hole's orientation changes following 
merger [Eq. (22)]. These results may be of particular observational interest. Jets launched by 
a black hole's spin should track its inclination: if the hole is kicked into a new orientation, the 
jet will "kink" . The angle of the kink should equal the hole's change in orientation. Applying 
our predicted dependence for the kink angle on the binary's parameters may provide some 
insight into the conditions of the hole's last merging, untangling a bit of its recent growth 
history. An abrupt change in inclination [such as discussed in Merritt & Ekers (2002)] 
requires a comparatively rare major merger. 

This work grew out of discussions at the KITP conference Black Holes: Theory Con- 
fronts Reality, Three Years Later, we thank the participants of that meeting for their stim- 
ulating input. We are also extremely grateful to Xinwu Cao for pointing out that the figure 
used in an earlier version of this manuscript contained important errors. SAH is supported 
by NSF Grant PHY-9907949; RDB is supported by NASA grants NAGW 5-2837, 5-7007. 

REFERENCES 

Bailyn, C. D., Jain, R. K., Coppi, R, & Orosz, J. A. 1998, ApJ, 499, 367 

Bardeen, J. M., Press, W. H., & Teukolsky, S. A. 1972, ApJ, 178, 347 

Begelman, B. C, Blandford, R. D., & Rees, M. J. 1980, Nature, 287, 307 

Binney, J. & Tremaine, S. 1987, Galactic Dynamics (Princeton University Press: Princeton), 
Chap. 3 

Blanchet, L. 2002, Phys. Rev. D, 65, 124009 



-9- 



Buonanno, A. & Damour, T. 1999, Phys. Rev. D, 59, 084006 
Colbert, E. J. M. & Mushotzky, R. F. 1999, ApJ, 519, 89 
Colbert, E. J. M. & Ptak, A. 2002, ApJS, submitted (astro-ph/0204002) 
Damour, T. 2001, Phys. Rev. D, 64, 124013 

Davis, M., Ruffini, R., Press, W. H., & Price, R. H. 1971, Phys. Rev. Lett., 27, 1466 
Elvis, M., Risaliti, C, & Zamorani, G. 2002, ApJ, 565, L75 

Ferrarcse, L. 2002, to appear in Current High-Energy Emission around Black Holes, Proc. 
2nd KIAS Astrophysics Workshop ed. C.-H. Lee (Singapore: World Scientific) (astro- 
ph/0203047) 

Gebhardt, K., Pryor, C, O'Connell, R. D., WiUiams, T. B., & Hesser, J. E. 2000, AJ, 119, 
1268 

Glampedakis, K., Hughes, S. A., & Kennefick, D. 2002, Phys. Rev. D 66, 064005 
Hughes, S. A. 2001, Phys. Rev. D, 64, 064004 
Kauffmann, G., & Haehnelt, M. 2000, MNRAS, 311, 576 

Kormendy, J. & Gebhardt, K. 2001, in AIP Conf. Proc. 586, 20th Texas Symposium on 
Relativistic Astrophysics, eds. J. C. Wheeler and H. Martel (Melville: AIP), 363 
(astro-ph/0105230) 

Lehner, L. 2001, Class. Quantum Grav. 2001, 18, R25 

Lifshitz, E. M. & Pitaevskii, L. P. 1980 Statistical Physics (Pergamon Press: New York), p. 
389 

Merritt, D. 2002, ApJ, 568, 998 

Merritt, D. & Ekers, R. D. 2002, Science, 10746881-0 (online at 
www.sciencemag.org/ cgi / content / abstract / 1074688vl) 

Miller, M. C. 2002, ApJ, submitted (astro-ph/0206404) 

Misner, C. W., Thorne, K. S., & Wheeler, J. A. 1993 Gravitation (Freeman: San Francisco), 
Chap. 33. 

Pathria, R. K. 1972, Statistical Mechanics (Pergamon Press: New York), p. 462 



-10- 



Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes 
(Cambridge University Press: New York) 

Rees, M. J. 1978, Nature, 275, 516 

Sasaki, M. & Nakamura, T. 1982, Prog. Theor. Phys. 67, 1788 
Small, T. A. & Blandford, R. D. 1992, MNRAS, 259, 725 

van der Marel, R. P. 2001, in Black Holes in Binaries and Galactic Nuclei. Proceedings of 
the ESO Workshop held at Garching, Germany, 6-8 September 1999, eds. L. Kaper, 
E. P. J. van den Heuvel, P. A. Woudt, 246 (astro-ph/0001342) 

Wilms, J., Reynolds, C. S., Begelman, M. C., Reeves, J., Molendi, S., Staubert, R., & 
Kendziorra, E. 2001, MNRAS, 328, L27 

Wilson, A. S. & Colbert, E. J. M. 1995, ApJ, 428, 62 

Yu, Q. & Tremaine, S. 2002, MNRAS, (in press) 



This preprint was prepared with the AAS IM^jX macros v5.0. 



-1 -0.5 0.5 1 -1 -0.5 0.5 1 

Cos i Cos i 



Fig. 1. — Spin a! of the remnant at mass ratio g = 0.1 (left) and q = 0.5 (center); right 
panel is a key. The axes are original spin a (vertical) and cosine of plunge inclination /i 
(horizontal). When q — 0.1, rapid rotation follows only if the original spin was large anc? the 
merger was nearly prograde (// ~ 1). About half of the parameter space yields a' < 0.4M. 
Higher mass ratio yields more rapidly spinning remnants: when q = 0.5, the plunge must be 
nearly retrograde for a' < 0.6, and no configuration yields a' < 0.5. 



